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Abstract 

The research herein studies the Langevin dynamics allowing for an exchange of energy between 
liquid crystals and the thermal environment. This dynamics leads to fluctuation and dissipation 
behaviors in the motions of liquid crystals, and therefore drives the system toward non-equilibrium 
evolutional processes. In particular, for the operations of liquid-crystal metastable systems, the 
fluctuation could allow an excitation (non-equilibrium) transition against energy barriers to the 
globally-stable state. Implemented with an actual case of liquid crystal tt configuration, this work 
statistically studies the non-equilibrium metastable transitions and shows the dependence of the 
transition-time on the correlations (of fluctuations). 
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I. INTRODUCTION 



In continuum models of liquid crystals jl-4|, which are based on phenomenological elastic 
constants and hydrodynamic transport coefficients, the equations of motion for the liquid 
crystal director are obtained by minimizing the free energy of the system. However, in 
liquid- crystal metastable systems [5], the energy-minimization process probably causes the 
liquid crystal state to be in a nearby local minimum of energy (weakly stable state) and pro- 
hibits the barrier transition to the observable global minimum state (strongly stable state). 
This process could result in a deviation between theoretical predictions and observations 
in the metastable systems. Accordingly, this paper considers the non-equilibrium Langevin 
dynamics js S] which allows an energy exchange between liquid crystals and the thermal 
environment, and introduces fluctuation £ (absorbing energy) and dissipation oc 7 (emitting 
energy) behaviors associated with the motions of liquid crystals. In principle, this is related 
to the fundamental fluctuation-dissipation relation (£(£i)£(£o)) oc j(ti ~ to), and exhibits 
noticeable influences in some liquid crystal systems. For instances, Barbero indicated that 
thermal fluctuations can lead to a re- normalization of the anchoring energy of nematic liq- 



uid crystals [9]. Galatola studied the Freedericksz transitions in nematic liquid crystals by 
thermal fluctuation approach n]]. Kelly showed the director fluctuations in nematic droplet 
with spherical harmonics method Ruhwandl presented the topological defects around 
a spherical colloid particle by Monte Carlo simulation 12j . Rehberg observed the spatial 
and temporal correlations of thermal fluctuations in liquid crystal systems [13|. Maclennan 
indicated a direct observation on the fluctuation of the liquid crystal director in his work 
[14 1 . In this work, we mainly introduce the influences of thermal fluctuations by Langevin 
dynamics to study the complete metastable transitions V — > (T, H) — > H in liquid- crystal 
7r configurations 15| . which still remains indefinite in real three-dimensional systems as we 
know. 

Based on Frank-Oseen elastic theory [lf|, we consider the spatiotemporal-correlated fluc- 
tuations £ as well as the corresponding dissipations oc 7 for the liquid crystal direc- 
tors in Langevin dynamics, since (I) the common non-correlated fluctuations £ n 17| con- 
ditioned by (£n(*i)£n(£o)) °c iksTditi — to) will cause a non-physical/divergent quan- 
tity in ultra-small time scale t\ — to ~ 0, while (II) the spatiotemporal-correlated fluc- 
tuations £ offers a more real algorithm to describe the experimental observations [3]. 
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Definitely, we consider the fluctuations £ related to an exponential correlation behavior 
(£(ri, ti)£(ro, to)) T^sTe - ' ^ °'e in the 1+1 (spatial+temporal) dimensions via 



the Ornstein-Uhlenbeck process 18| , and thereby characterize the fluctuation behavior by 



correlation time r and length A as experimental observations [13J. Implemented with an 
actual case of the metastable systems in the liquid-crystal 7r configuration 15], this pa- 
per fulfills a statistical numerical analysis over an ensemble of systems with 1000 identical 
iterations, and elucidates the non-equilibrium metastable transitions. 

This article is organized as follows. Section 2 reviews Frank-Ossen elastic free energy 



using Q tensor representation 16|, and structures the fundamental equation of motion for 
liquid crystal directors. Section 3 further introduces the Ornstein-Uhlenbeck process [18( for 
the Langevin dynamics to consider an actual spatiotemporal-correlated fluctuations £ (r, t) . 
An update formula of directors for the numerical calculation is then identified. Section 4 
numerically analyzes the influences of fluctuations and dissipations in the liquid crystal it 
configuration, and depicts the associated non-equilibrium metastable transitions, i.e., about 
how a initial V state transit through the locally-stable T and finally into the globally-stable 
H state in the real three-dimensional systems. The main conclusions are then discussed 
in section 5. Two more experiments as well as the corresponding theoretical works are 
appendixed and show the influences of fluctuations and dissipations even in usual cells of 
Liquid-Crystal Displays (LCDs). 

II. MODELING LIQUID CRYSTAL DIRECTOR 

In elastic continuum theory, the local aligned sample may be regarded as a single liquid 
crystal, in which the molecules are on the average aligned along the direction defined by a 
director n. The Frank-Oseen strain free energy density using tensor representation for liquid 



crystal directors can then be described by 



I q(2) ! q(2) 

fs = ^(K 33 - K n + 3K 22 )-^ + ~{K n - K 22 )-±- 

-i ^(3) r (2) 

+ l (K 33 -K 11 )^ r + qo K 22 ^ (1) 
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with 



G[ •* — Qjk,lQjk,l C4 ■* — £jklQjmQkm,l (2) 
^2 ^ = Qjk,kQjl,l Gq ^ = QjkQlm,jQlm,k (3) 



in which 



Qjk = S (rijUk - ^jkj , (4) 

Qjk,i = -^jr»- i>M e ( 5 ) 

Here, i^n, i^ 2 2, and K 33 are the splay, twist, and bend elastic constants, respectively. q = 
27r/po associates the pitch p of chiral liquid crystals, e.g. Cholesterics. Term Q is the 
tensor order parameter, and n is the liquid crystal director with unit length. S is the scalar 
order parameter, and indicates the ideal isotropic and nematic phases by the values S = 
and S = 1, respectively. Ejki is the Levi-Civita symbol. 5jk represent the Kronecker delta. 
Moreover, the convention of summing over repeated indices is used herein. The electric free 
energy density can be expressed in terms of the Q tensor as 

f e = + \e A eV d V^, (6) 

with 

e = g , Ae = e||-e_i_, (7) 

dV 

Vj = j e {x,y,z}. (8) 

Here, e is the LC dielectric tensor with elements e± and e\\ for dielectric coefficients perpen- 
dicular and parallel to the LC director. 

To solve equilibrium states of liquid crystal directors with applied potentials, the system 
free energy, which is obtained by integrating the Gibbs free energy density f g — f s — f e over 
the volume, is minimized. Accordingly, the Euler-Lagrange equations for electric potentials 
and director components under condition \n\ = 1 can thus be obtained 



= - \fg\v = d i [( £ ± S ij + ( £ \\ - e±)n i n j )d j V] , 
= - [fglm + Xn i ie {>, 



(9) 
(10) 



Here, A is the Lagrange multiplier for constraint \n\ = 1. By applying the chain rule 
[fglm = if^Q k^T = [/»]q-* ' S (njdki + nkdji), we obtain the alternative definitions in 
tensor representation for Equation (TTOT) 

[f 9 ]n x = ' 2 "^ f:, o,, ■ :l ■ /</o, 

+n z s([f g } Qxz + [f g ] Qzx ), (11) 

[/ 9 ]n y = 2n v S [fg\ Qyv + n * S (lf 9 ] Qxy + \fg\gj) 

([/Jq,. + Wo J > (12) 



+%5 ([/ 9 ] Qyz + [/ 9 ] Q J (13) 

Solving Equations (J2]) and ( TTOT) simultaneously, the potential distribution and the director 
configuration for equilibrium status can now be evaluated. It is noted that, in general, \fg\q fc 
is not equal to [f g }q k . for cases of j ^ k. For the analyses of dynamic time evolutions on the 
liquid- crystal directors, however, we consider the time-dependent /non-equilibrium Langevin 
equation. In the Langevin dynamics which includes the stochastic torque (fluctuation) aQ 
as well as the dissipation torque k for directors, the general rotational Langevin equation 
under a external torque F can be read as [17| : 

d 

n—rii = -Fi + a£, (14) 



Here the constant a is related to k by Einstein's relation a = a/2/^bT, and ( is the normal- 
ized white noise with dimension s -1 / 2 . For the liquid-crystal system, to simply evaluate the 
dissipative torque for the liquid crystal director, the volume of the liquid crystal director is 
treated as an effective sphere with diameter d, and thereby the director will suffer a dissi- 
pative torque k = n^/d 3 by Stokes relation. An arranged Langevin equation corresponding 
to Equation (fT0|) is then obtained: 

7^ = " If,}* + ^ 

= - La, + 1 (15) 

Details for this part will be discussed in the next section. Moreover, the equation of motion 
for common use [l6[ can be estimated by considering the statisticly steady director (n^), i.e. 
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ignore the (temperature-dependent) thermal fluctuations £ 



9\{m) 



(16) 



with an effective order parameter < S < 1 in [f g ]/ n .y Note that the function of the 
Lagrange multiplier A has been replaced here by re- normalizing the director |n| = 1 at each 
time step [16]. 



III. SPATIOTEMPORAL CORRELATED FLUCTUATION FUNCTIONS 

Based on the above-mentioned Langevin dynamics for liquid-crystal systems as the Equa- 
tion (115p . we now detail the spatiotemporal-correlated stochastic torque (fluctuation) £ as 
well as the corresponding dissipations for the liquid-crystal directors. Mainly we consider 
the Ornstein-Uhlenbeck process in 1+1 dimensions [is| to generate the spatiotemporal- 
correlated stochastic torque £(z,t) with: 

«0M)> = (17) 

<«*■«)«*■*■» = ^C?""^ < i8 > 

Here, X z and r are the characteristic correlation lengths in the spatial (z) and temporal 
(t) dimensions, respectively, ks is the Boltzmann constant with a value of 1.38 x 1CT 23 
NmK -1 , d is the effective size of the liquid crystal directors, and T is the temperature of 
thermal environment. The Equation ([TBI) can be further understood by its power spectrum 
P(z ,oo) = i+t 2 lu 2 wm ch results from the temporal Fourier transform of the correla- 
tion function (£(zo,t)£(zo,0)). It is indicated that a longer correlation time r suppresses 
the higher-frequency (higher-energy) energy exchanges and leads to a smaller fluctuation 
amplitude [ ^fix rl 1- Likewise for the spatial Fourier transform of io)£(0, t )), i.e. 
P{k,t ) — XtfC \ + \Zk 2 ' a l° n g er correlation length A z implies the smaller- momentum (longer- 
wavelength) energy exchanges, and causes a longer-ranged but smaller fluctuation strength 

\ 7k B T 1 ll 
l2nd?T A J • 

To generate this real-space fluctuation function £(z,t) to evaluate the motion of liquid 
crystals in Equation (IT5|) . an exact recursion relation with temporal spacing A t and spatial 



spacing A 2 in the finite difference method 



18} is applied: 
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£(z,t) = -exp(-^--!-)z(z-A s ,t-A t ) 




x 1 - exp ( ) Ar] (z, t) 



(19) 



Here, r](z,t) is spatiotemporal Gaussian white noise with (r}(z,t)rj(z' ,t')) = S ZjZ /8 tt t', and A 



spatiotemporal parameters and the characteristic of white noise function (r](z,t)) = in 
Equation ( fl9l) . £(z, t) is independent of its initial value and evidently meets the condition in 
Equation (IT7|) . For the condition of fluctuation correlations (£(z, t)£(zo, to)) in Equation ffl8|) . 
numerical statistics averaged over 10 3 iterations (each iteration includes the data number 
iV = 100 in spatial/temporal space) are fulfilled to verify the validity of Equation (fT9j) . The 
numerical results by Equation ( fl9l are shown to compare well with the analytical results 
of Equation ( fl8l) in Figures 1(b) and 2(b) for temporal and spatial arguments, respectively. 
Figure 1(a) shows £(zo,t) of real temporal space in one of the iterations and characterizes 
a r-period fluctuation function. Figure 2(a) then presents to) of real spatial space with 
coherent length X z . 

Implemented with Equation (JT9~|) for the real-space fluctuation function t), an updated 
formula by discretizing Equation ( TT5l) can be obtained: 



for the equation of directors in 1+1 dimension systems. This algorithm includes the in- 
fluences of thermal baths and supplies a theoretical foundation for non-equilibrium liquid 
crystal dynamics. 



is the fluctuation amplitude [ 2 JdVA J^ - ^ ue ^° ^ ne ex P onent ially decaying correlation for 



rii(z,t + At) - rii(z,t) 



At 



(-[/<?U + £M))> i^{x,y,z} 



(20) 
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FIG. 1: Analyses of numerical fluctuation function £ for fixed spatial state: z = z$ and = 2A Z . 
(a) Fluctuation function in real temporal space for t = At (thin red line) and r = 5Aj (bold blue 
line), (b) Temporal correlation function for r = A( (symbol □) and r = 5A^ (symbol O)- Solid 
lines indicate the corresponding analytical results. 




50 100 5 10 

(a) z (b) z-z 



FIG. 2: Analyses of numerical fluctuation function £ for fixed temporal state: t = to and r = 2Af. 
(a) Fluctuation function in real spatial space for X z = 2A Z (thin red line) and X z = 8A Z (bold 
blue line), (b) Spatial correlation function for X z = 2A Z (symbol □) and X z = 8A Z (symbol Q). 
Solid lines indicate the corresponding analytical results. 
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FIG. 3: Liquid crystal configurations and the corresponding Gibbs free energy in ir configuration: 
bend (V), twist (T), and homogenous (H) equilibrium states. 
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IV. NUMERICAL ANALYSES OF NON-EQUILIBRIUM DYNAMICS 



To concretely understand the influence of the thermal fluctuations on the liquid-crystal 



metastable systems, we now consider an actual case with tt configuration [15j. As shown in 
Figure 3, with the Frank-Oseen elastic theory, bend (V), twist (T), and homogenous (H) 
states form a metastable status. A transition barrier exists between the T and H states, 
leading to a local minimum localization on the T state for the time evolution from the 
V state. (Here the red, green, and blue arrows indicate the spatial x, y, and z directions, 
respectively.) In this work, we mainly explain how the initial V state can transit through the 
locally-stable T state, and finally into the globally-stable H in the long run as experimental 
observations. To statistically evaluate the influences of the fluctuation and dissipation, we 
consider an ensemble (10 3 iterations) of liquid crystal cells for each set-up conditions of 
parameters. The associated parameters are defined as: 7 = 50mPa ■ sec, A" n = 13.6pN, 
K 22 = lO.OpN, K 33 = U.7pN, £|| = 10, e± = 5, and T = 300K. The thickness of the cell 
is set as 4.5 \im. No chiral agent is added. The pre-tilt angle between the liquid crystal 
directors and the alignment (green) layers is set to 5°. For the finite difference algorithm, the 
spatial grid is A z = 0.15/xm and the temporal grid is A t = 0.01ms. Besides, the effective 
director size d is treated as about 10 1 — 10 2 nm for this study, which is 1 — 2 order larger than 
the liquid crystal molecules. For each iteration of the ensemble, the V state is initialized 
by applying a bias 15V, which lasts for a long enough period (200ms in our cases), and 
the voltage is turned off at t > in order to observe the evolution of the liquid crystals. 
Figure 4(a) illustrates the transition rate T H = l/(tv^n) versus fluctuation amplitude A, 
where (£y_>#) represents the average time of V — >■ H transitions of the ensemble. In the 
Figure 4(a), in the absence of the thermal bath (at A = OPa), the process of minimizing- 
energy as the Equation f|T6|) drives the dynamics of the liquid crystals, and consequently the 
initial V states in all of the 1000 iterations are found to follow the free-energy-decreasing 
trace into the locally-stable T state as the blue solid-line segment form V to T in the 
Figure 4(b). No iteration of transition from V to H is found in our simulation in this 
case. This results in the infinite transition time = 00 for V to H state, or the zero 

transition rate Yh = 1/(V-Kff) = as in Figure 4(a). (Note that Figure 4(a) shows the 
log-scale plot for clear appearance, and the infinite time is realized by a cur-off time 10 5 ms 
in real numerical calculations.) As the fluctuation amplitude A increases (A 7^ OPa), the 
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liquid crystals contain not only the elastic free energy but also the finite fluctuation energy; 
Or alternatively, the liquid crystals suffer an extra stochastic torque £ in addition to the 
fundamental elastic torques. This allows the initial V state in some/few iterations of the 
1000 iterations have the possibility to deviate from the free-energy-decreasing trace, and are 
found to proceed another non-equilibrium path (as the green dashed curve in Figure 4(b)) 
to reach the final H states, although the paths of the most iterations still remain on the free- 
energy-decreasing trace (red dashed-dotted curve in Figure 4(b)) to be localized at T states 
in our calculations. This causes the transition time (ty^u) of some of the 1000 iterations 
to be finite and thereby the finite statistic transition rate T H = 1/ {ty^u) as in Figure 4(a). 
Consequently, it is found that larger fluctuation amplitude A can cause more iterations of 
the ensemble to proceed the non-equilibrium path (V — > H; green dashed curve in Figure 
4(b)) and raise the transition rate Th = l/ifv-tii) as in Figure 4(a). Here, Figure 4(b) 
illustrates quantitatively the energy evolution traces of the transitions V — > T and V — > H 
for A = 100 Pa, r = A t , and \ z = A z . For the large fluctuation amplitude (A — > oo 
Pa), however, the overwhelming fluctuation torque leads to strong- fluctuated liquid crystal 
orientations, and thereby the liquid crystal configurations (V, T, and H) cannot be well 
identified. This is far away from phenomenological observations and is beyond the scope of 
this work. In addition to these analyses of the transition dynamics on the aspect of time 
and energy as in Figures 4(a) and 4(b), respectively, We further analyze the corresponding 
time evolutions of the directors to help to understand the transition dynamics. 

Figure 5 purposes to briefly characterize the liquid-crystal profiles of the above-mentioned 
transitions first, and the detailed time-evolution traces of the corresponding liquid-crystal 
directors are then depicted in Figure 6. The characteristic director profiles of the transition 
V — > T, which corresponds to the red dashed-dotted curve on Figure 4(b), is characterized on 
Figure 5(a). The characteristic director profiles of the transition V — > H, which corresponds 
to the green dashed curve on Figure 4(b), are found to have two possible time-evolution 
traces in our simulations and are characterized by So and Si in Figures 5(b) and 5(c), 
respectively. In Figure 6(a) (corresponding to Figure 5(a)) the orientations of the individual 
liquid-crystal directors draws the bend (V) state at t — and denoted by (n x , n y ) — >(0, 0), i.e. 
\n z \ — > sj\ — ri% — riy = 1. By the time going, these directors then proceed (along the red 
dashed-dotted energy trace in Figure 4(b)) to the locally-stable T state, and consequently 
describe the transition V — > T over the long run. Here, an alternative anti-clockwise-twist 
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FIG. 4: (a) Transition rate of V — > H as a function of fluctuation amplitude A with characteristic 
correlation scale r = At, X z = A z (mark □), r = At, X z = 5A Z (mark +) , and r = 5At, X z = A z 
(mark <>)■ In which the horizontal axis corresponds to the d = 93nm, 47nm, and 31ram for A = 50, 
100, and 150-Pa, respectively for r = At, X z = Az, and T = 300-fT. And for the case of r = At, 
X z = 5Az and r = 5Ai, X z = Az, A = 50, 100, and 150-Pa corresponds to d = Alnm, 20nm, 
and lAnm, respectively, (b) Energy evolution traces of alternative non-equilibrium transitions for 
V -> T (dashed-dotted curve) and V -»■ H (dashed curve) for A = 100 Pa (T = 300 if with d = 50 
ram), r = A 4 , and A 2 = A z . 

profiles is analogous to the presented clockwise-twist one, and is ignored in the context. 
Similar anti-clockwise/clockwise configurations occurred on the transition V — > H is ignored 
below as well. Figures 6(b) and 6(c) numerically show two possible time-evolution traces 
of directors (n x ,n y ) corresponding to the transitions V — > So — > H and V — > Si — > H, 
respectively. Note these two transitions (V — > So — > H and V — > Si — > H) show an 
indistinguishable energy trace as the green dashed curve in Figure 4(b) in our simulations. 
All these construct the fundamental understandings to help to study the complete V — >■ H 
transition in real three-dimensional liquid-crystal n systems as below. 

In the following section, we numerically investigate an extension study of the three- 
dimensional liquid- crystal metastable system to give a result corresponding to the actual 
experimental observations [151 ] . It is emphasized that the spatial/temporal grid size (A z 
and A 4 ) is herein assumed to be much larger than the coherent length/time (X z and r), 
such that the fluctuations for each time step at the individual site are almost independent, 
i.e. £ (z,t) = Ar](z,t) in Equation (jl~9l) . This thereby allows a straightforward extension 
to the three-dimensional condition £ (x, y, z, t) = A'rj (x, y, z, t) for the studied system. The 
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FIG. 5: Characteristics of liquid crystal configurations of alternative non-equilibrium transitions 
for (a) V — > T process, (b) V — > Sq — > H process, (c) V — > S± — > H process viewed from +x (left) 
and — y (right). 




(a) n 10 (b) n 1 ' (c) n 1 



FIG. 6: Evolution traces of liquid crystal director (n x ,n y ) corresponding to (a) V — > T process, 
(b) V — > So — > H process, (c) V — > Si — > H process. In which the directors at z = d/9 (mark •), 
z = 2d/9 (mark +), z = 3d/9 (mark x), z = Ad/9 (mark o), z = 5d/9 (mark □), z = 6d/9 (mark 
0), z = 7d/9 (mark <), z = 8d/9 (mark [>) are shown. 

relevant parameters of the 7r configuration are similar to the above-mentioned, except for 
A z = 0.45 fim as well as the considerations of optical transmittances [l9j related to practical 
observations. The associated optical parameters are set as: the incident wavelength is 0.55 
fim, ordinary refractive index of liquid crystals is n = 1.5, extraordinary refractive index 
of liquid crystals is n e = 1.58, and the absorption axes of the upper and lower polarizers are 
located at 4> = 45°, and 135° in the x — y plane, respectively. The numerical results for the 
metastable transition V — >■ T,H without (A' = 0) and with (A' = 100) thermal couplings 
are shown in Figures [7] and El respectively. In the case of A' = in Figure El the initial 
V state follows the energy-minimizing process in the elastic continuum theory and leads to 
the localization at the T state as described above. However, the observable H state exhibits 
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(a) V state (b) T state (c) forbidden H state 

att=0ms at t= 1000 ms 




(d) LC profile 
at t^O ms 



(e) LC profile 
at t= 1000 ms 



(f) LC profile 
of H state 



FIG. 7: Optical transmittances of (a) the initial V state at t = ms, (b) the stable T state at 
t = 1000 ms, and (c) the forbidden H state for the V —> T,H transition without thermal coupling 
A' = 0. The corresponding liquid crystal directors (LC profiles) at z = d/2 are showed as below 
plots (d), (e), and (f), respectively. 

significant differences of liquid crystal configurations to the T state in Figures [7(e) and [7(f) 
and is forbidden by the energy barrier (Figure 11(b) ). In the case of A' = 100 in Figure [5J 
thermal fluctuations as well as the associated correlations offer the possibility of the non- 
equilibrium V — > H transitions and cause the T + H transient states in 100 ms after the bias 
is dropped to zero (Figure [8(b)). Thereafter, the elastic dynamics describing the disclination 
motions [15( drives a final H stable state as in Figure [8(c). This gives an explanation for 

the overall V — > H process. Overall, with thermal fluctuations, a straightforward analysis 

I I 

corresponding to the observed state transitions of n- configuration in Reference [15( can be 
realized as in Figure (jHJ). 

V. CONCLUSIONS 

In this work, based on the elastic continuum theory, we have introduced the Langevin 
algorithm with the Ornstein-Uhlenbeck process to realize the non-equilibrium dynamics of 
liquid crystal systems coupled weakly to the thermal environment. In which the Ornstein- 
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(4 LC profile <e) LC profile (f) LC profile 

at Mi rns att=100ms attlOOOms 



FIG. 8: Optical transmittances of (a) the initial V state at t = ms, (b) the transient T + H 
state at t = 100 ms, and (c) the stable H state at t = 1000 ms for the V —> T,H transition with 
thermal coupling A' = 100. The corresponding liquid crystal directors (LC profiles) at z = d/2 are 
showed as below plots (d), (e), and (f), respectively. 




FIG. 9: Transition processes of 7r-configuration under the bias sequence. In the field-on period, the 
homogenous (H) state transforms to the stochastic symmetric (Hs) and asymmetric (Ha) states, 
and further proceed to the bend (V) state under higher biases. In the field-off period, the bend 
(V) state first relaxes to the twist (T) and homogenous (H) states as described in the context, and 
recovers to the overall homogenous (H) configuration in the long run. 
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Uhlenbeck process supplies a connection to the natural spatial and temporal correlation 
scales of thermal fluctuations. All these allow one to simulate a non-equilibrium time evolu- 
tion and give physical insights into the processes taking place in the studied systems. It is 
noted that another thermal hopping ( hop over the potential barrier, i.e. a direct transition 
T — > H in Figure 11(b) ) dynamics can be set up into this system through strong coupling 
cases, but go beyond this study's main goal. Definitely, this work provides one possible 
explanation for the complete transition V — > H in the tt configurations and supplies the 
corresponding real-time simulation results for the experimental observations 151 ]. 



Appendix A: Fluctuation-induced Vortex/anit vortex in VA cells 

In this appendix, we briefly show the fluctuation- induced (stochastic) vortex/antivortex 
through experimental and theoretical works, and the associated fluctuation amplitude A' can 
be estimated as well. The used structure of the VA (vertical alignment) liquid crystal cells 
is depicted in Figure (ITUj) . The associated LC parameters are defined as: 7 = 135mPa ■ sec, 
if 11 = 13.6piV, K 22 = lO.OpN, K 33 = U.7pN, e\\ = 3.6, and e± = 7.3. Here, the fluctuation 
amplitude A' are decided by fitting the density and period of vortexes/antivortexes with 
experimental observations, and obtained as A' = 80 — 125Pa. These amplitude value corre- 
spond to the fluctuated angle of the director around 5° for studied cases in our simulations. 
The pre-tilt angle between the liquid crystal directors and the alignment (xy) plane is set to 
90° as the experimental setup. For optical analyses, the incident light wavelength A is set 
as 0.55 /j,m, n = 1.5, n e = 1.6, and the absorption axes of the upper and lower polarizers 
are located at = 0°, and 90° in the x — y plane, respectively. Besides, a long enough 
pre-simulation (lasting 100 ms) is fulfilled to simulate the natural/initial LC status for the 
sequential numerical simulation. Numerical results as well as experimental observations for 
optical transmittances are illustrated in Figure [TT] by time. It indicated that the thermal 
fluctuations can introduce stochastic vortex/antivortex as observed in experiments HIT a), 
while without fluctuations a regular time evolution of LC profiles is predicted as in Figure 
[TTTc). Definitely, these allow a further analysis on the emergence and annihilation behav- 
iors of vortex/antivortex associated the spatial and temporal correlations of fluctuations 
while beyond the scope of this work. 
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FIG. 10: Structures of the VA liquid crystal cells with patterned ITO (indium tin oxide) electrodes. 
The cell volume is set to be 50um x 50um x 4um in width x length x height (x x y x z) dimensions. 
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{3) Experimental Observations 



mm 




(b) Simulation with thermal fluctuations 





(c) Simulation without thermal fluctuations 



FIG. 11: Optical transmittances of the VA cell under bias 5V in (a) experimental observations, (b) 
numerical analyses with fluctuations, and (c) numerical analyses without fluctuations. In which 
(a) and (b) show the comparable emergence and annihilation behaviors of vortex/antivortex by 
time, (c) shows the regular time evolutions of LC profiles on transmittances. 

Appendix B: Effects of Fluctuations on LC response time 

In this appendix, not only the fluctuation-induced vortex/antivortex as described at Ap- 
pendix A, the thermal fluctuations are shown to have appreciable influences on the LC 
response behaviors (against biases). In the meantime, this influence supplies an alternative 



16 



method to quantitatively evaluate the fluctuation amplitude A'. The used structure of the 
MVA (multi-domain vertical alignment) liquid crystal cells is depicted in Figure ffl2|) . in 
which the upper ITO (indium tin oxide) has fine-slit patterns and the lower ITO stacks with 
protrusion layers for the tilted-alignment purposes. The associated parameters are referred 
to these at Appendix A; while the fluctuation amplitude A' are decided to be 115 Pa by fit- 
ting the response time as below. Numerical results as well as the experimental observations 
for the optical transmittance at t=16 ms and V=5 V are shown in Figure [131 Normalized 
transmittances as a function of time at biases V=4, 5, and 6 V are then depicted by numeri- 
cal simulations with fluctuations in Figure [HJa), numerical simulations without fluctuations 
Figure [T4T b). and experimental observations in Figures (H^a) and (b) for comparison. Com- 
pared to the common understandings: larger bias/field causes faster LC response in field-on 
process as shown in Figure [T47 b) (solid curves; without fluctuations), thermal fluctuations 
as well as experimental observations give another story herein. Besides fast LC response in 
field-on process, large bias/field cause simultaneously dense vortex spots as shown in Figure 
[131(a). These observations agree with the numerical results if the fluctuations are included 
as Figure fTST b) . Definitely, the sequential reliefs of the entangled vortexes could cost more 
time to approach the final stable status, and thereby could reduce the liquid crystal re- 
sponse time. This can be observed in Figure [T4T a): compared between these curves with 
V=4 and 5 V, the response time with larger bias (5V) is indeed faster than that with lower 
bias (4V); however, compared between these curves with V=5 and 6V, the response time 
with larger bias (6V) is slower than that with lower bias (5V). In Figure fT4Tb). however, the 
numerical simulations without fluctuations could not help to understand these experimental 
observations. Substantially, this gives an example to indicate the significance of fluctuations 
in some real liquid crystal systems. 
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Lower [TO 
with Protrusions 



FIG. 12: Structures of the MVA liquid crystal cells, in which the upper ITO has fine-slit patterns 
and the lower ITO stacks with protrusion layers. The cell volume is set to be 160um x 200um x 4um 
in width x length x height (x X y X z) dimensions. 




(a) Experimental (b) Simulation with (c) Simulation without 
observation thermal fluctuations thermal f I uctuations 



FIG. 13: Optical transmittances of the MVA cell with bias 5 V at t=16 ms by (a) experimental 
observations, (b) numerical simulations with fluctuations, and (c) numerical simulations without 
fluctuations. In which (a) and (b) show the stochastic vortex/antivortex spots, (c) shows the 
regular transmittance pattern. 




FIG. 14: Normalized transmittances as a function of time at biases V=4,5,6 V for (a) numerical 
simulations with fluctuations, and (b) numerical simulations without fluctuations. The experimen- 
tal results are added to both analyses for comparison. 
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